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ABSTRACT 

We present observations of the starburst galaxy NGC 7714 with the Infrared Spectrograph IRS 1 
on board the Spitzer Space Telescope. The spectra yield a wealth of ionic and molecular features 
that allow a detailed characterization of its properties. NGC 7714 has an HII region-like spectrum 
with strong PAH emission features. We find no evidence for an obscured active galactic nucleus, 
and with [Nelll] / [Ne II] ~ 0.73, NGC 7714 lies near the upper end of normal-metallicity starburst 
galaxies. With very little slicate absorption and a temperature of the hottest dust component of 
340 K, NGC 7714 is the perfect template for a young, unobscured starburst. 
Subject headings: dust, extinction, galaxies: individual (NGC 7714, Arp 284), galaxies: starburst 



1. INTRODUCTION 

NGC 7714 is a peculiar barred-spiral galaxy at 45° in- 
clination. For cz — 2798 km/s the galaxy is at distance 2 
of 39.6 Mpc . Together with its post-starburst compan- 
ion NGC 7715 it forms the interacting system Arp 284, 
which is the result of a recent (100 — 200 Myr), off-center 
collission between the two disk galaxies l|Struck fc S mith 
12003ft . With its compac t, UV-bright nuc l eus, N GC 7714 
has been classified by iWeedman et alJ l|1981fl as the 
proto-typical starburst galaxy. The central region of 
about 330 pc has been the site of active star formation at 
a rate of about 1M@ yr _1 for some 10 8 years. However, a 
recent significant increase in the star forma tion rate made 
it the dominating source of the UV flux l)Lancon et alJ 
2001). 

With its strong Hell A4686 line 

l|Oonzalez-Delgado et al.l TT995T . NGC 7714 has been 
classified as a Wolf-Rayet galaxy. In fact, the optical 
and UV spectra indicate a population of about 2000 
Wolf-Rayet and 20000 O-type stars, suggesting a 
fairly young age of the pre sent starburst of 4 — 5 Myr 
IjGarcia- Vargas et al.l 119971). This is i n agre ement 
with earlier st u dies by IWeedman et al.1 IJ1981T and 
iTaniguchi et all (|1988[) who estimated 10 4 05-type 
stars from the Balmer and Bracket-7 line fluxes, and 
masses of ionized gas of 3.0 x 10 6 M Q , and 1.9 x 10 6 M Q , 
respectively. The radio-continuum and the X-ray 
luminosity of 6 x 10 40 ergss -1 require 10 4 supernova 
remnants in a volume of radius 280 pc IjWeedman et al.1 
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1981). ITaniguchi e t al. (1988) found definite evidence 
for a star burst-driven bipolar winds from the nucleus, 
nearly perpendicular to the disk plane. 

Although most of the activity is concentrated in the 
nucleus, NGC 7 714 as a whole is experiencin g intense 
star formation ( Gonzalcz-D elgado et al.lll995l) . Mid-IR 
imaging with ISOCAM IjO'Halloran et al]l2004ft revealed 
a strong source at the nucleus surrounded by slightly ex- 
tended emission out to abo ut 30" (5.7 kpc) i n diame- 
ter. From the IRAS fluxes l)Surace et al.ll2004|) and the 
above distance we calculate a total infrared luminosity 
Of Lg-lOOOAim = 5.6 x 10 10 L Q . 

The strong optical emission lines, H/3 and [O III ] , 
show no signs of b road e mission l)Weedman_e j_aL|jl98lL 
ITaniguchi etTal] l|1988ft L Recently. iSoria fc M otch 
(2004) found two compact X-ray sources with XMM- 
Newton. One of them coincides with the starburst nu- 
cleus, has an X-ray luminosity of Lx ~ 10 41 ergs -1 , 
and shows the spectrum of a thin thermal plasma with 
a power-law (point-source) contribution. The variability 
in the power-law component hints at the presence of ei- 
ther a hidden low luminosity AGN or an ultraluminous 
X-ray source l|Soria fc Motcbll2004[) . 

In this letter we describe new observations of 
NGC 7714 with the Spitzer Space Telescope. These ob- 
servations are part of an IRS guaranteed time program 
to obtain high signal-to-noise spectra of a large sample 
of different classes of nearby galaxies, that can be used 
for comparison with more distant systems (Devost et al. 
2004, Higdon et al. 2004a, Houck et al. 2004a). 

2. OBSERVATIONS AND DATA REDUCTION 

We observed NGC 7714 with the Infrared Spectrograph 
(IRS) (iHouck et alJl2004bft on board the Spitzer Space 
Telescope l|Werner et al.l 120041) . The observations are 
part of the IRS guaranteed time program. The data 
were taken during the first IRS campaign in nominal 
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operations on 17 December 2003 using the standard 
IRS "Staring Mode" Astronomical Observing Template 
(AOT). The dataset consists of observations with all 
four IRS modules: 2 cycles x 14s in "Short-low" (SL, 
AA = 5.2 - 14.5/an, R ~ 64 - 128), 2 cycles x 14s in 
"Long-low" (LL, AA = 14.0 - 38.0/mi, R ~ 64 - 128), 4 
cycles x 30s in "Short-high" (SH, AA = 9.9 - 19.6/zm, 
i? - 600), and 2 cycles x 60s in "Long-high" (LH, 
AA = 18.7-37.2/mi, R ~ 600). Each cycle yields two ex- 
posures at different nod positions along the slit. The slits 
were positioned relative to the reference star SAO 128273 
using an "IRS high accuracy peak-up" in the blue filter 
band. 

The data have been pre-processed by the Spitzer Sci- 
ence Center (SSC) data reduction pipeline version 9.5 
(Spitzer Observer's Manual, chapter 7 3 ). The two- 
dimensional "basic calibrated data" (BCD) constituted 
the basis for further processing. The BCD frames have 
been individually inspected by eye, and hot pixels, as 
well as very negative pixels, that have not already been 
flagged by the pipeline, were masked manually. The for- 
mer could easily be identified by flipping between im- 
ages at different nod positions. Tests have shown that 
masked pixels introduce artificially high noise in the spec- 
tra, and hence we interpolated isolated masked pixels by 
their nearest neighbor pixel values along the dispersion 
direction. Next we averaged the frames from the same 
nod positions (in the case of SH where four frames per 
nod position were available we have used the median in- 
stead). Since the low-resolution, long-slit modules con- 
tain two sub-slits, each integration provides a "free" sky 
spectrum - mainly zodiacal emission - in the adjacent 
subslit. We computed the median sky and subtracted it 
from the low-resolution frames. 

Further processing was done within the IRS data 
reduction and analy sis package SMART, version v. 4 
( Hiedon et al. 2004lJ - a powerful IDL package for spec- 
tral extraction and spectral analysis. The high-resolution 
spectra were extracted with "full aperture" extraction 
along the diffraction orders. The ends of each orders 
where the noise increases significantly were manually 
clipped. Within any one module the individual orders 
matched remarkably well and required no further fine- 
tuning, with the exception of SH order 11 which was 
scaled up by 5%. Finally the two nod positions were av- 
eraged. The low-resolution spectra were extracted using 
SMART'S "interactive column extraction" , which is sim- 
ilar to the method used in the SSC pipeline. The 3rd 
"bonus" order in SL and LL has not been included. 

At the current stage of calibration, there remains a sig- 
nificant mismatch between some modules in both low- 
and high-resolution spectra. This mismatch is most 
likely d ue to either (i) extended source emission (as seen 
e.g. by O'Ha lloran et al.1 (J2£)04)), or (ii) pointing errors 
that lead to flux losses that are most significant at short 
wavelengths (narrow slits). We also note that the low- 
resolution slits are narrower than the slits of the high- 
resolution modules. At this early stage of the Spitzer 
mission we have no means to identify the real cause. 
However, the LH and LL spectra agree very well with 
each other, and a comparison betw een the IRAS 25/im 
flux of 3.15 Jy <|Suraceet alJ 12004) and the LL flux in 

3 http://ssc.spitzer.caltech.edu/documcnts/som/ 



the IRAS filterband of 2.61 Jy yields good agreement, 
taking into account that the IRAS "aperture" is much 
bigger and is even more susceptible to extended emission. 
Adopting the LH flux densities as correct, we scaled SH 
up by 36% to achieve an excellent match of the contin- 
uum fluxes around 19/im. We note that the flux in the 
[S III] line, which lies in the overlap region between SH 
and LH, is then the same in both modules. Similarly, the 
low-resolution spectra from different modules had to be 
scaled to match. The LL 2nd order was scaled up by 17% 
to match the LL 1st order; the SL 2nd order was scaled 
up by 40% to match the SL 1st order; and finally the re- 
sulting SL spectrum was scaled up by 17% to match the 
LL spectrum. As a result of this approach the low- and 
high resolution spectra agree very well with each other, 
and - even more important - we determine from the low- 
resolution spect rum in the IRAS 12u m band a total flux 
of 0.47 Jy, while iSurace et alJ l)2004|) found 0.56 Jy from 
IRAS - again good agreement! Figures Q] and El show the 
final IRS high- and low-resolution spectra, respectively. 
The properties of the ionic and molecular features were 
obtained from single Gaussian fits to the high- and low- 
resolution spectra within SMART; the results are listed 
in tables ^ and El 

3. RESULTS 

3.1. The starburst properties of NGC7714 

Based on optic al and near-IR spectroscopy 
IjLancon et alJ 1200 lj) . the IS OCAM filter band ra - 
tios and PAH emission features (O'Halloran et al. 2004), 
it has long been argued that the nucleus of NGC 7714 
is a site of intense starburst activity. However, given 
its low extinction, the unusually high infrared lumi- 
nosity and recent rep orts of a possible, hidden AGN 
IjSoria fc Motchl I2004J) more precise diagnostics are 
needed to understand the processes in the nuclear 
region. The high signal-to-noise mid-IR IRS spectra 
provide the ideal diagnostic tools to characterize the 
nature of the underlying power source in more detail. 

Table ^ lists the properties of the strongest fine- 
structure lines in the high- resolution spectrum (fig. ^). 
We find 10.8 x lQ- 20 Wcm- 2 for the [Nell] line while 
iPhillips et all IJ198J) found (7 ± 2) x 10- 20 Wc m - 2 from 
the ground. For the [S IV] line we dete rmine 1.8 x 
10- 20 Wcm- 2 while lO'HaHoran et al.1 ll20f)l miote (3.2± 
0.8) x 10~ 20 Wcm- 2 based on ISOPHOT-S low-resolution 
spectra. The S(0), S(l), and S(2) lines of molecular hy- 
drogen are marginally detected and will be discussed in 
a subsequent paper. 

Of particular interest is the [Ne III] / [Ne II] ratio as 
it is a good measure of the hardness of the interstel- 
lar radiation field, which is mainly determined by the 
most massive stars. We find [NeIII]/[Ne II] =0.73 ± 
0.05, corresp onding to an effectiv e temperature of about 
3.8 x 10 4 K l|Oiveon et all 12001. Based o n the "ISO- 
SWS Starburst Sample", iThornlev et~all <|2000|) stud- 
ied the [Ne III] / [Ne II] ratio for 27 galaxies, with val- 
ues ranging from 0.05 (M83) to 12.0 (IIZw40), and a 
median value of 0.26. The starbursts of the nearby 
ISO sample with [Ne III] / [Ne II] ratios closest to our 
value are NGC 3690 (0.71) and the interaction region in 
NGC 4038/39 (0.84). The origin of these variations is 
still controversial (Rigby & Rieke 2004, Thornley et al. 
2000). Even within the same galaxy the [Ne III] / [Nell] 
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Fig. 1 — IRS high-resolution spectrum of NGC 7714. The total 
integration times are 240 seconds for both SH and LH (see discus- 
sion in the main text). 
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Fig. 2. — IRS low-resolution spectrum of NGC 7714 (see discus- 
sion in the main text). The total integration time is 56 seconds per 
subslit. 
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EMISSION FEATURES IN 
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ID 




A obs 

M 


E pa 




Flux b 


S/N a 


EW a 


-Vest 


A obs 




Flux 


EW 


Relative 


Mod a 




H 


[eV] 


[io- 


- 20 Wcm~ 2 ] 






H 


H 


H 


[10 


- 19 Wcm" 2 ] 


[ium] 


strength 15 




[S IV] 


10.51 


10.61 


34.8 




1.79 ±0.04 




92 


0.03 


6.2 


6.31 




4.60 ± 0.44 


0.50 


0.49 


lo 


[Nell] 


12.81 


12.94 


21.6 




10.82 ±0.48 




135 


0.11 


7.7 


7.80 




9.42 ± 0.39 


0.70 


1.00 


lo 


[Nellll 


15.56 


15.63 


41.0 




7.95 ±0.11 




111 


0.11 


8.6 


8.71 




1.74 ±0.28 


0.17 


0.18 


lo 




SIII] 


18.71 


18.89 


23.3 




9.01 ±0.10 




131 


0.08 


11.2 


11.33 




2.21 ±0.28 


0.17 


0.23 


hi 




SIII] 
Sill] 


33.48 


33.77 


23.3 




12.61 ±0.51 




37 


0.11 


12.7 


12.81 




1.66 ±0.16 


0.14 


0.18 


hi 




34.82 


35.13 


8.15 




10.45 ±0.36 




48 


0.09 


16.4 


16.55 




0.33 ±0.03 


0.03 


0.04 


lo 



















Note. — Emission line properties obtained from a single Gaussian 
fit to the high-resolution data. 

a EP — Excitation potential, EW — Equivalent width (observed), S/N 
— signal-to-noise ratio. 

The uncertainties quoted for the line fluxes throughout this letter are 
the errors from the line fit and do not include the calibration uncertain- 
ties. 



ratio can var y significantly as rec ently seen with the IRS 
in NGC 253 Ipevost et all [2004). We will address this 
important issue in a subsequent paper once we have 
observed a larger sample of starburst galaxies. How- 
ever, with its large number of very luminous Wolf-Rayet 
stars, it is not surprising that the [Nelll] / [Nell] ratio in 
NGC 7714 lies above the value found in most starburst 
galaxies. 

Although [OIV] has an excitation potential of 54.9 eV, 
faint emis sion has b een seen in many starburst galax- 
ies IjLutz et al.lll998|) and recently with the IRS in th e 
ultra-luminous IR-galaxy UGC 5101 i|Armus et all2 004). 
We also find a marginal detection of the [O IV] line with 
(9.9 ± 7.2) x 10" 21 Wcm" 2 . Since the location of the 
line coincides with a noisy pixel, we take this value 
as an upper limit. Hence, with [OIV] /[Nell] < 0.09 
and a relative PAH(7.7/^m) /continuum(7.7/im) w 1.30, 
NGC 7714 falls below the regio n occupied by ULIRGs 
in fig. 5 of iGenzel et alJ 1^)98). Furthermore, the ab- 
sence of the [NeV] line - we place a 3c upper limit of 
3.9 x 10 -21 Wcm - emphasizes the pure starburst na- 
ture of NGC 7714. On the basis of the infrared emission 
lines, therefore, we find no evidence for an obscured ion- 
ized region associated with an AGN, which implies that 



Note. — PAH band emission strength, For the measurement of 
the 12.7/im feature we first removed the [Nell] line. 

a lo — fit to low-resolution data, hi — fit to high-resolution data 
Relative to the 7.7/im feature. 



the va riable X-ray source reported by ISoria fc Motchl 
( 2004) is either from strong shocks associated with very 
recent supernova activity, or emission from a HMXB. 
There is growing evidence that both possibilities are 
intimately associated with massive young star clusters 
IjGaoet al]l2003JK 

3.2. The dominant PAH features 

Below 13/jrn, the spectra a re dominated by st rong 
emission from PAHs (see e.g. iPeeters et alJ l|2003j) for 
a recent overview), with the strongest at 6.2/im, 7.7/im, 
8.6/im, 11.2/im, 12.7/im, and 16.4/xm listed in tabled 
More features (e.g. at 5.7/mi) appear to be present 
but will not be discussed here. Some of those fea- 
tures have also been seen in the ISOPHOT-S spectra 
which cov ered the 5 — 12/xm ran ge at a resolution of 
R = 90 IjO'Halloran et alJ I2004J) . The total flux we 
determine for the strongest feature at 7.7/im agrees 
with the ISOPHOT-S measurement of (8.84 ± 1.10) x 
10- 19 Wcm- 2 . 

Shortward of 13^m the spectrum is qualitatively very 
similar to the well-know starburst galaxy M82. How- 
ever, unlike M82, which has very strong PAH emission 
with respect to the warm continuum from very small 
grains and plausibly some extinction due to silicates at 
9.7/im IjForster-Schreiber et alJl2003|) . the spectrum of 
NGC 7714 shows a rather smooth continuum and no ev- 
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idence for silicate absorption. 

With the high signal-to-noise provided by Spitzcr-IRS 
it now becomes possible to expand the comprehensive 
study of fainter PAH features seen in Galactic sources 
to a large sample of extragalactic sources. While the 
strength of e.g. the 6.2/mi PAH / continuum ratio is re- 
markably constant over a wide luminosity range in Galac- 
tic sources, others, e.g. the 11.2/im feature, which is 
linked to neutral PAHs, get quickly reduc ed in stronger 
i onizing radia t ion fields ( Honv et alJl200l|) . 

Honv et al. (2001) show that the 12.7/xm band corre- 
lates well with the CC stretch mode at 6.2/zm for Galac- 
tic sources ranging from YSOs, HII regions, and reflec- 
tion nebulae to planetary nebulae (their fig. 5). With 
log 10 (/ 6 .2//ii.2) = 0.32 and log 10 (Ii2.r/Iii.2) = -0.12, 
NGC 7714 falls right on their correlation line. We con- 
clude that the physical conditions in NGC 7714 are sim- 
ilar to Galactic HII regions. In fact, a c omparison 
with the ISO-SWS spectrum of Ml 7 l|Peeters et al.l2004|) 
yields a perfect overall match longward of 15/xm. 

3.3. Dust and Extinction in NGC 7714 

One of the most remarkable attributes in Figure is 
the lack of silicate absorption in NGC 7714. This im- 
plies that the observed spectral shape of the mid-IR con- 
tinuum is defined purely by the continuum emission of 
the hottest dust directly heated by the young, massive 
stars. If so then NGC 7714 might be a well-needed tem- 
plate for fitting relatively unobscured starbursts at low 
and high redshifts. Is NGC 7714 the perfect, unobscured 
starburst? 

From t he hydrogen recombinat ion lines Pa/3, Br7, 
and Bra, iPuxlev fc Brand I IJ1994J) determined Ay = 
1.8 ± 0.7 mag for a point source model with obscur- 
ing screen, and Ay = 3.9 ± 1.7 mag if sources and 
dust are homogeneously mixed. For an extinction law 
of Ag,Q^ m /Ay = 0.1 and if Ay = 3.9 we find that 
79.6/jm = 0.36 and hence the continuum around 9.6/im 
would be suppressed by a factor of at most 30%. In 
the case of an obscurin g screen, which was favored by 
IPuxle v & Brand (1994), the suppression is only 15%, 
which is in good agreement with the shape of our spec- 
trum. 

The long wavelength baseline available in our spec- 
trum of NGC 7714 makes possible a very accurate de- 
termination of the continuum, which allows sophisti- 
cated fits for the dust emission which is responsible for 
this continuum. Figure shows a three component fit 
to the low-resolution spectrum. Two components are 
warm and cool astronomical silicate/graphite (Si/Gr) 



dust with an MRN distributio n llMathis et al.l | 1977[) and 
the optical cons t ants g iven bv lDraine fc Lee 1 1) 19841) and 
lLaor fc Draine I 1^)93). The third component is astro- 
nomical PAHs for which we us e the absorption coefficient 
of Sicbenmor gen et alJ 1)20011) . We also allowed for the 
possibility of an intervening, absorbing Si/Gr dust screen 
but found that r 9 .7 Mm < 0.2 (consistent with the above 
results). The Si/Gr dust components have temperatures 
of approximately 153 K and 60 K, and the PAH compo- 
nent has a temperature of about 340 K. PAH emission 
clearly dominates below 9/xm while the Si/Gr mix fills 
in the "PAH trough" at 10.5/im. The Si/Gr component 
then dominates at longer wavelengths. A warm and a 
hot component are clearly present in NGC 7714. Due 



IRS low-resolution spectrum 
cool Si/Gr component (T-60K) 

n Si/Gr component (T=1 53K) 
hot PAH component (T=340K) 




Observed wavelength [|im] 



Fig. 3. — IRS low-resolution spectrum of NGC 7714 with the 
various fit components overlaid (see discussion in the text). Strong 
emission lines were excluded from the fit. 



to the absence of silicate absorption, we have confidence 
that the continuum is directly observed, without being 
obscured by cooler, intervening dust. This means that 
we have achieved a direct measurement of the hottest 
dust in a prototypical galaxy. This determination is im- 
portant because the dust temperature is often used as a 
discriminant between starbursts and AGN. 

This work is based [in part] on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under NASA contract 1407. Support for this 
work was provided by NASA through Contract Number 
1257184 issued by JPL/Caltech. 
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